function [alpha] = ode(omega, theta)

R = 100;
H = 10;
n = 4;
g = 9.807;

alpha = ((n.^3) * (H.^2) * (omega.^2) * sin(n * theta) - n * g * H) * cos(n * theta) / ((R.^2) + (n.^2) * (H.^2) * (cos(n * theta))^2);

end

